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We compute, in the large iV limit, the quark potential for TV = 4 supersymmetric 
SU(N) Yang-Mills theory broken to SU(Ni) x SUfflz). At short distances the quarks see 
only the unbroken gauge symmetry and have an attractive potential that falls off as 1/L. 
At longer distances the interquark interaction is sensitive to the symmetry breaking, and 
other QCD states appear. These states correspond to combinations of the quark-antiquark 
pair with some number of VK-particles. If there is one or more VF-particles then this state is 
unstable because of the coulomb interaction between the VF-particles and between the W's 
and the quarks. As L is decreased the VF-particles delocalize and these coulomb branches 
merge onto a branch with a linear potential. The quarks on this branch see the unbroken 
gauge group, but the flux tube is unstable to the production of VF-particles. 
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1. Introduction 

The conjectured link [Q between large iV M = 4 supersymmetric Yang-Mills and ten 
dimensional supergravity on AdS^ x S5 has led to a number of interesting predictions. 
In particular one has been able to compute the large iV behavior for the attractive po- 
tential between heavy quarks and antiquarks, both at zero temperature and finite 
temperature 0,^,0, and for a quark and a monopole 0, collections of dyons, monopoles 
and quarks and flux tubes 0. In these calculations, one looks for static string con- 
figurations that correspond to a separation L and potential energy E between the heavy 
particles. 

For a quark-antiquark pair at zero temperature there is a unique string configuration 
for a given separation L 0|| . At finite temperature, one finds that there are some values of 
L where two solutions exist . The branch with the lower energy is the coulomb branch, 
where at short distances, the quark and anti-quark experience an attractive coulombic 
force. The other branch is a "confining" branch, in that the force is linear in L, at least 
for small L. The physical meaning of this branch is unclear. 

There are other systems where more than one branch exists. As we will see below, an 
SU(N) gauge theory broken to Yli SU(Ni) with gNi large will have many branches. In fact 
as L diverges, the number of branches diverges. In this paper we will study the breaking 
of SU(N) to SU(N/2) x SU(N/2). In this situation some of the string configurations can 
be solved by quadratures, or exactly in various limits, and the qualitative features can be 
generalized to other breakings. 

The physics of these extra branches can be deduced by looking at the energy density 
along the static string configuration, and also by taking the quark separation to infinity. 
In particular we find static string configurations that correspond to states consisting of 
a quark, an antiquark and n massive W particles, with the latter transforming in the 
{N/2, N/2). The energy density along the string concentrates appropriately: there are two 
(infinitely) massive quarks at the ends and n equal and evenly spaced mass concentrations 
in between. In the limit of infinite separation the total energy of the system becomes the 
sum of the BPS quark masses plus the BPS masses of n VF-particles, while the energy 
density reduces to that of the quarks and n evenly spaced 5-functions for the VF-particles. 
The static string configurations in this limit also decompose into a concatenation of such 
configurations for a quark, an antiquark and n intermediate W's. As the quark separation 
for the n VF-particle configuration is decreased the branch becomes unstable to small 
perturbations because of the attractive coulomb potential between the W particles and the 
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quarks. As L is decreased further the VF-particles delocalize and the unstable branches 
merge onto a branch with a potential that is linear for large enough L. 

Turning this around, for large enough quark separation there is an unstable flux tube 
configuration that has a linear potential between the quarks. This is an unbroken gauge 
configuration in that the quark and anti-quark are both in the fundamental of SU(N). 
The flux tube can then destabilize by popping out VF-particles that are then subject to 
the coulomb instability. Eventually, a stable configuration is reached where we have a quark 
and an anti-quark in (N/2, 1) or (1, N/2) representations of the broken group. Finally, as 
L is decreased even further, these stable branches merge onto the unbroken gauge branch. 
For L below this critical value this is the only string solution. 

In section 2 we discuss some general features for static string configurations and discuss 
how the study of such configurations reduces to the study of particle motion in a potential, 
V. We use Maupertuis' principle to argue that the force between the quarks is always 
attractive and has a simple relation to the classical particle's energy. We then specialize 
to geodesies that correspond to an SU(N) gauge theory broken to SU(N/2) x SU(N/2), 
for which the potential, V, is bipolar and symmetric. In section 3 we consider string 
trajectories that lie on the symmetry axis. We show that there is a branch of trajectories 
that runs up a saddle and which is coulombic for small separations and confining for large. 
We also show that there are infinitely many branches that merge onto this confining branch. 
All but the lowest energy branches are unstable to small perturbations and eventually 
become coulombic as the separation is increased. In section 4 we consider trajectories that 
are off the symmetry axis. In this case, the quark and anti-quark have unequal masses. 
The branches split into groups of three, except for the lowest energy branch which is 
isolated. The higher energy branches have a minimum quark separation. Pushing the 
quarks inside of this separation leads to a finite release of energy. Section 5 contains some 
final comments. 

2. Metrics, Geodesies and the Quark Potential 

2.1. Generalities 

The general form of the metric for a configuration of parallel, static D3-branes in ten 
dimensions is: 

ds 2 = H' 1 / 2 dz^dz v n^ + H 1 / 2 dy m dy m , (2.1) 
where /j, = 0, . . . , 3 and m = 4, . . . , 9, and H is a function of the y m alone. 
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Our purpose is to use the ideas first mooted in [JJ, and further deveioped in 0- 
II , to study the quark-antiquark potential. This means that we consider static string 
configurations that begin and end on a single brane at a very large (and ultimately infinite) 
distance, y = y m ax, from the rest of the branes in the configuration defined by H. Finding 
such static string configurations is equivalent to finding geodesies in the spatial 9-metric: 

ds 2 = Hiy)- 1 dz a dz a + dy m dy m , (2.2) 

where z a , a = 1,2,3, runs over the spatial coordinates of the D3-brane. Because of the 
symmetry in the z a , we can choose coordinates so that the geodesic runs in the z 1 = z 
direction, with z 2 = z 3 = 0. Let A be an affine parameter for geodesies in the metric ( |2.2| ), 
then the translational symmetry implies that 4f = const H(y). We will normalize the 
affine parameter by taking the constant equal to one, and hence 

Tx = H ^ < 2 ' 3 > 

Using this conservation law in the geodesic equations for y m (X) yields: 

^ + 29^=°- (2 ' 4) 

We are thus studying the motion of a particle in the (repulsive) potential V = \H. We 
will denote the classical energy of this particle by: 

M = U'^Y + \H(V) ■ (2.5) 



2 V d\ J 2 
The energy of the quark-antiquark pair is given by 



and the three-dimensional separation, L, of the quark-antiquark pair is given by: 

H(y(X)) dX . (2.7) 



The energy density along the string can be readily deduced from ( plf) and (|2.7|) , with the 
density p given by 

p = w m ■ (2 - 8) 
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We ultimately wish to study E as a function of L, and seek out coulombic and confining 
behavior. One can gain some very useful insight into E(L) by applying some of the variants 
of variational principles associated with particle motion, most particularly, Maupertuis' 
principle [ID|. Consider a general variation of the action S(-S^q) about a solution of the 



Euler-Lagrange equations, but for which we let the endpoints and the time for the path 
vary: 

5S = (pfSq f -H f St f ) - (pi-Sft-HiSti) , (2.9) 

where the subscripts i and / denote initial and final quantities, while p is the particle 
momentum, and H is the Hamiltonian. We will consider variations for which p ■ Sq = 0, 
and in partcular variations in which the endpoints do not move. Since the system has time 
translation invariance with a conserved energy, [i, we have 5S = —fiSt. For such a system, 
it immediately follows that the particle path is obtained from the variational principle 
5Sq = 0, where So is the abbreviated action: 

S Q = S + n{t f -ti) = Jp-dq. (2.10) 

In doing this variation one considers all paths with a fixed energy, fj,, and one must express 
the integrand of So in terms of li, q and dq only. For example, for a particle moving in a 
potential V, one has: 



S = J v / 2(/j-V)dq-dq. (2.11) 

Consider now, the variation of the path with respect to the energy, \x. The time taken 
along the path will change in response, and so we have from ( |2.1U| ): 

5S = 5S + (t f -U)Sn + fiSt = (t f -U)6n, (2.12) 

and so <ga = (tf-U). 

Returning to the original problem, we want to consider how the quark energy, E and 
separation, L, change as we vary the geodesic, and in particular, vary the classical particle 
energy, fx. From ( |2.6|) , (|2.7|) and (|2.12|) , one has 



8E = -= ( a SX + IXSa) 

rra' \J2jji 



(2.13) 



5L = 25 / V = 5(S - 25) = XSjj + 2jj8X . 



From this one immediately sees that 



8E = 1 —= SL , or ^ = 1 —= . (2.14) 

2ixa'^ dL 27ra'V 2 7* 
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Thus the force between a quark and anti-quark is determined by the energy of the classical 
particle. 



There are several important consequences of this. First, the force is always attractive: 
something that is obvious from QCD, but not at all obvious from the string perspective. 
In particular, if one considers E and L as a function of some other parameter, as one 
does in the thermal case |5|,[| , one sees that E and L must reach a minimum or maximum 
at the same point. Hence if one of these functions "turns over" then so must the other. 
Such simultaneous turning points of E and L are crucial to the cusp-like behavior E(L) 
at bound state transitions. This will be discussed in more detail in section 4. 



On a more fundamental level, to compute the quark potential one considers a family 
of geodesies going out to the brane at infinity. A confining family is thus characterized 
by a family of such geodesies for which the classical particle energy is constant. One can 
directly verify this in a special case: the situation where the family consists of a particle 
that is approaching a local maximum, or saddle point from below, and turning around 
and reversing its course just before the maximum, or saddle. Suppose the potential at 
the maximum or saddle is po, and consider a particle coming in with energy po — e. For 
the particle to approach the maximum, the time taken diverges as t ~ —log(e), and so 
the inter-quark energy E, and length L, also both diverge as —log(e), which means that 
E ~ L. 



It is also interesting to consider the behavior of E and L for particles that approach, 
or are released from a point where the potential is diverging as V(q) ~ q~ r . One can 
easily see that L ~ p~*^ and E ~ p~^ as one approaches the singularityS. One thus finds 

2 

E ~ L~' F ^ , which is consistent with ( |2.14|) , and shows Coulombic behavior only for r = 4. 



We can also find the energy density along z, where we see that the density is p ~ p~ x ^ 2 
for the part of the string near the peak and p ~ p 1 ^ 2 for the part that is away from the 
peak. Hence, very energetic particles correspond to string configurations where the energy 
density is highly clumped. 



1 One must, of course, subtract off the quark mass to get this fall off in the energy. 
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2.2. Multi- centered D3-branes 



In QXJ] a single group of iV D3-branes was considered, and the function H was given 



by: 



H 



1 + 



4ngNa 



12 



= v / y m y m . 



(2.15) 



In the large iV limit, one then has the anti-de Sitter metric 



ds* 



U 2 dU 2 

— dz^dz u rj„ u + R 2 — + R 2 dn 2 , 



(2.16) 



where U = r/a' , and R = (AngN) 1 / 4 . It was also pointed out in |l]] that one could equally 
well consider separated groups of D3-branes, and that for large N and large U these would 
also limit to anti-de Sitter space. 



The multi-centered D3-brane solution is then given by taking |Tl| , |12 

4irgN a a' 2 



H(y) = 1 + 



\y - Vc 



;2.i7) 



where N a is the number of branes located at the point y a . As in [Q , we consider the limit 
of the metric in which a' — >• with r/a' finite, and in which all of the iV a are large (but 



N a /N(, is finite). This means that we may drop the 1 in (|2.17| ). The IIB superstring in this 
background should then describe a Higgs phase of an iV = 4 super symmetric SU (N) gauge 
theory, with iV = ^2 a N a , where the gauge group has been broken to Yl a SU(N a ). The 
Higgs vevs are represented by the separations of the groups of D3-branes. In particular, 
at large r = ^/y rn y m , the function H in (|2.17| ) takes the form fl2.15|) with iV = N a , and 
thus at large distances, and for large N a , the multi- centered solution limits to the anti-de 
Sitter solution corresponding to SU(N). 

Our major focus will be the breaking of SU(N) to SU{M) x SU(M), M = N/2, and 
upon geodesies in the corresponding double-centered metric. By rotating the coordinates, 
we can take the potential V 



\H to be: 



12 



V(x,y) = — 



R 4 



+ 



R 4 



((x - a'4>) 2 + y 2 ) 2 ((x + a'4>) 2 + y 2 ) 2 



;2.18) 



where we have set y 1 



y. To make contact with the conventions of other authors 



in whose work a' has been scaled away and in which the variables x and y have 
dimensions of mass, one makes the rescaling x — > a'x, y — > ct'y and A — > a' 2 A. 
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We will analyze the geodesic structure extensively in the next section, but here we 
note some general features. First, for large r = \/ x 1 + y 2 , one has V = o/2.R 4 /(2r 4 ), and 
provided the geodesic stays out at large values of r, (r >> a'(j)) the properties of such a 
geodesic will not differ significantly from those discussed in [§||§. In the language of gauge 
theory this means that the symmetry breaking mass scale of SU(N) to SU{M) x SU(M), 
which is proportional to (p, is much less than the energy scale of the quark-antiquark 
interaction {1/L) 1 and so the interaction behaves as it would in the unbroken SU(N) 
phase. 

At the other extreme, there are two geodesies that represent single quarks: That is, 
the geodesies start on the single D3-brane and terminate on one or other of the groups of 
M D3-branes. They traverse no distance in the z-direction (the constant in the conserved 
quantity is zero) and run with infinite classical particle speed. Such geodesies represent 
quarks in the (M, 1) or (1,M) of SU(M) x SU(M) and the quark mass is simply the 
distance in y-space from y ma x to the location of the relevant group of branes. 

There is a third important geodesic: it also does not move in the z-direction, but runs 
directly from one group of M D3-branes to the other. The corresponding classical particle 
runs from one peak to the other through the interconnecting saddle. This can be thought of 
as representing the massive VF-particle coming from the broken gauge symmetry, and lying 
in the (M, M) or (M, M) representation depending upon the orientation of the geodesic. 
The mass of this W is given by ( |2.6|) , but without the potential term (the constant of 
integration for the z-motion is zero) and is thus given by mjy = (fi/n. 

Returning to the geodesies that correspond to a quark-antiquark pair, we will see in 
the next section, that as we separate the pair, i.e. increase L, the geodesic drops closer 
to the double- centered core. Initially the pair is in an unbroken phase, but at a critical 
separation there is a phase transition and one can begin to distinguish the (1, M) and the 
(M, 1). Associated with this transition, one also finds geodesies that correspond to states 
with the quark, anti-quark and several (M, M) or (M, M) VF-particles. 

3. Analysis of Trajectories 

Using the potential derived in the preceeding section, we can now analyze the quark- 
antiquark potential. We first consider trajectories that come in at, or infinitessimally close 
to an angle of 90 degrees with respect to the SU (M) x SU (M) scalar axis. The advantage 
of doing it this way is that the H-i symmetry is not broken by the U(l) scalar expectation 
value. Moreover, the lowest trajectory can be solved by quadratures. 
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3.1. The geodesic on the symmetry axis 



Consider the trajectory that is aimed exactly at the saddle of the potential. The 
trajectory then stays at 90 degrees since there is no lateral force on the particle. From 
the potential in ( |2.18|) , we see that we have an effective potential in y only, while the x 
coordinate stays fixed. Using ( |2.7| ) and (|2.6D , the length L on the D3 brane between the 
quarks is 

L = 2 f d\V(0,y(X)) , (3.1) 



(3.2) 



while the quark energy is given by the integral 

E = h I dX V(dxy) 2 + 2V(0,y(X)) 



where the term inside the square root is twice the conserved energy for the particles motion. 
(We have passed to the conventions of by recaling x — > a'x, y — > a'y and A — > a' 2 \ as 
described earlier.) 

Using (|2.3|) and conservation of the particle energy, we see that as a function of z, y 
satisfies 

4V 2 (0 : y)(d z y) 2 = C-2V(0,y) , (3.3) 

where C is an integration constant. Hence, for a particle tractory that starts at y = oo 
goes down to y = yo and then turns around and goes back to y = yo, we have 



V(0,y) 



y/2V(0,y )-2V(0,y) 



dy 



2R 2 (<j> 2 +y 2 ] 



dy 



(3.4) 



(y 2 + <P 2 )V(y 2 + <P 2 ) 2 -(yo 2 + <P 2 ) 2 



The above integral is a combination of elliptic integrals of the first and third kind. In 
terms of generalized hypergeometric functions the length between the quarks is 



(2n) 3 / 2 R 2 



(r(i/4))2v^+^ 



3-^2 



13 3.15. 

4' 4' 4' 2' 4' 



(yo 2 + (P 2 ) 2 



<t>\T(l/4)Y 



AM 2 {(f) 2 +J/Q 1 



, ■ x x , p i 3 5 5. 3 7. y 

^ ' ■■" ' ,2^| 3 2 I 4' 4' 4' 2' 4' ty Q 2 _|_ ^A2^|2 



(3.5) 



The quark energy in (|3.2|) is found using ( p73|) and (|3.3| ). Hence, we find that 

"OO 



E 



1 

7T 



y 2 + <^ 2 



(3.6) 
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This integral is divergent because the energy includes the quark bare masses, which diverges 
as the position of the brane is taken to infinity. We can find the finite potential piece by 
cutting off the integral at y ma x and subtracting off y m ax/n~, the sum of the quark masses. 
Alternatively, we can regularize the integral in ( |3.6| ) by replacing the numerator of the 
integrand by (</> 2 + y 2 ) s and analytically continuing to s = 1. In either case, the result is 
the same and the quark potential energy is 



E 



^vV 2 + yo 2 

(T(l/4))2 



13.1 

4' 4' 2 



(yo 2 + (p 2 ) 2 
r(i/4))V 



16n 2 (<p 2 + ?/o ; 



Ft 



I 5.3. 

4' 4' 2' („,_2 



(yo 2 + 



\2\2 



The asymptotic behavior of L for large and small yo is 

(2n) 3 / 2 R 2 



L 
L 



(r(i/4)) 2 

V2R 2 



yo- 1 + O(l/y 3 ] 



log(yo) + 0(1) 



yo » 
yo « <P , 



while the asymptotic behavior for the potential is 



(r(i/4)) 2 

log(y 



yo + O(l/y ; 
+ 0(1) 



2/o » 



yo « 



7T 



Hence we see that for small L the quark potential has the coulombic behavior 

~ (r(l/4))4 
while for large L the quark potential has the confining behavior 



E 



1 



! L = 



L . 



(3.7) 



(3. 



(3.9) 



(3.10) 



(3.11) 
This is 



27rv/4ipT r 2 R 2 

The coulombic behavior is the same as that found for the unbroken phase 
not surprising; at short distances the energy scale is much larger than the W mass, so one 
expects the potential to approach the coulomb potential for the unbroken case. 

However, for large L one would still expect to find a coulombic potential. This is 
because when the energy scale drops well below that of the W mass, one would expect the 
VF-particles to become irrelevant, and so one should find the coulombic behavior for an 
SU(M) (M = N/2) super Yang-Mills theory. Hence there should be other geodesies that 
describe this behavior. 
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3.2. Geodesies around the symmetry axis 



In order to find such geodesies, we start by looking for geodesies that are very close 
to the straight trajectory. Indeed, to lowest order in perturbation theory we can assume 
that the motion in the y direction is the same as for the straight trajectory, and take the 
motion in the x direction to be infinitesimally small. 

The equation of motion in the x-direction is then is approximately 



d 2 x d 2 V(x,y) 



d\ 2 dx 2 



_ 2R\H 2 - (y(A)) 2 ) 

(02 + (y(A))2) 4 V-U) 

x=U 



where y(X) corresponds to the motion of a straight trajectory. Hence if y(X) 2 > 5</> 2 then 
the motion in the x direction will accelerate exponentially away from the symmetry axis, 
x = 0. However, if y(X) 2 < 5</> 2 then the there is a restoring force, and the motion is 
oscillatory. Consider the geodesies that start at x = and y = y ma x — * oo and return 
to the same point. Such geodesies cannot exist unless they go close enough to the origin 
for the x-motion to experience a restoring force (i.e. have y(X) 2 < 5</> 2 ). It follows that 
if L « R 2 /(p then the only geodesic that starts and finishes on the symmetry axis is the 
one that runs down the symmetry axis and back. There is thus a phase of this theory 
in which the Z2 symmetry is unbroken, and the quark-antiquark pair feel only the full 
SU (N) gauge theory. 

Consider a geodesic whose initial velocity is such that x is still small when it enters 
the region with y(X) 2 < 5<fi 2 , and begins to oscillate. For the geodesic to begin and end 
on the symmetry axis, x = 0, the number of half-period oscillations has to be an integer. 
Hence, these modes will be quantized. We first examine the situation where there are a 
large number of oscillations, which means that many of the oscillations are very close to 
the saddle. We want to find the difference in energy between the n and n + 1 half-periods. 
The motions are almost identical, except for an extra half-oscillation near the saddle. For 
this part of the motion y can be approximated as a constant y ~ 0, and so using ( |3.12| ) we 
find that the period for this extra half oscillation is 

AA = JL 4> 3 . (3.13) 



Hence, using ( |2.7[ ) , ( |2.18| ) and ( |3.11| ) , the increase in quark distance between the appearance 



of n and n + 1 oscillations is 



AL = AXR 4 (f)~ 4 = ^2= (f>~ 1 ■ (3.14) 

v/To 
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and the change in the potential energy is 



AE = 



7T <f> 



(3.15) 



2VlO 7T 



Thus, a trajectory with an extra half-oscillation appears when the energy is increased by 
an amount that is slightly less than half the W mass, mw = 4*/ 71 - 

Once we have isolated a trajectory with a fixed number, n, of small half-oscillations, 
we can increase L and find the asymptotic behavior for this trajectory. As L increases, 
the trajectory will move further and further up the two peaks, or closer to the singularities 
of V, until it reaches the limit where the trajectory is the concatenation of a straight line 
coming in from (x, y) = (0, —00) to the peak at ((f), 0) and then n — 1 lines going back 
and forth through the saddle between the peaks and finally a straight line from one of the 
peaks back to asymptotic infinity. Figures 1 and 2 show the L behavior for geodesies with 
2 and 6 half-oscillations. Hence the total energy of this configuration is the sum of the 
quark masses plus n — 1 times the W mass 0/7T. In other words, as L — > 00, we find that 
we are left with the two heavy quarks plus n — 1 VF-particles. 




Fig. 1: Geodesies with two half oscillations overlayed on a contour map of the poten- 
tial. As L is increased, the geodesic approaches the concatenation of the BPS geodesies 
of two quarks and a W. The potential peaks are at the upper corners. The straightline 
trajectory in the figure is the small oscillation limit for two half-oscillations. 
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Fig. 2: Geodesies with six half oscillations. As L is increased, the geodesic approaches 
the concatenation of the BPS geodesies of two quarks and a five Ws. 



It is also instructive to reverse this argument, and start with the heavy quarks sepa- 
rated by a large distance but with n Ws in between. As L becomes finite, the trajectory 
moves off the peaks and has n + 1 smooth half-oscillations. Hence this corresponds to a 
situation with heavy quarks and n intermediate VK-particles. As L is decreased the oscilla- 
tions become smaller and smaller, until we hit a critical distance L n , where the oscillations 
disappear all together and the trajectory becomes our original straight-line trajectory. As 
we decrease L further, the straightline trajectory becomes shorter until it coincides with 
the trajectory with n oscillations at L = L n _i. Obviously, this trajectory corresponds to 
the case with heavy quarks and n — 1 VF-particles. Therefore, when L = L n -i, the binding 
energy for the last W among n VF-particles is the W mass <f>/7t. 

We can make a reasonable estimate of the strength of the coulombic potential on the 
branch with n VF-particles. If L is large, we can see from the figures that the oscillations 
get very close to the peaks. Moreover, each oscillation reaches roughly the same distance 
e away from the peak with its velocity becoming quite small. Hence we can use the results 
of ]2],[| to estimate L and E. For each half-oscillation there is an ascent and a descent of 
a peak. Hence, L function of e is approximately 



where there is a factor of n + 1 coming from the n + 1 half-oscillations and a factor of 
l/y/2, since the effective gauge group at each peak is SU(N/2). Likewise, the potential 
energy is 




(3.16) 




(3.17) 
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and hence E as a function of L is 



E 



~, — ; — ', — m -Lj 



(T(l/4))< 



(3.18) 



This coulombic behavior can be understood as follows: we can assume that there are 
n intermediate VF-particles equally spaced between the quarks. In the large N limit each 
W feels a coulombic attraction only with its two neighbors (the VF-particles on the end 
feel an attraction to the neighboring quark), where the strength of this interaction is the 
coulombic quark anti-quark interaction for gauge group SU(N/2) but for a separation of 
L/(n + 1). Since there are n + 1 interacting pairs, we expect the total coulomb potential in 
( |3.18| ) . Since the VF-particles can move, we should expect that these branches are unstable 
because of a coulomb instability. We will discuss this further in the next subsection. 

Figure 3 shows the branch of straight geodesies, which becomes linear for large L, 
and the first three coulomb branches. The derivatives dE/dL are continuous where the 
coulomb branches attach to the saddle branch, but d 2 E/dL 2 is not. The lowest coulomb 
branch is the quark-antiquark branch with no intermediate VF-particles. This branch is 
doubly degenerate, since a trajectory on this branch goes up one of the two peaks. The 
second coulomb branch is a quark and antiquark with an intermediate VF-particle. There 
is only one branch since these trajectories are TLi symmetric. Finally the third coulomb 
branch, which is also doubly degenerate, has two VF-particles between the quark-antiquark 
pair. 
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Fig. 3: Plot of E vs. L for the saddle branch and the first three coulomb branches. 
E is plotted in units of the W mass, mw, and L is in units of R 2 /mw- The first and 
third branches are doubly degenerate. The second branch is a TL^ symmetric branch. 
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It is instructive to compare the linear behavior of the straight trajectories to the 
coulombic behavior of ( |3.18|) when L = L n and n is large. In the linear regime ( |3.11| ) 
yields 

f L _ w „ L671 -! . ( , 19) 

where m\y = <f>/ir. For large n (|3.14|) implies L n »s ^=^ _1 , and hence the energy on the 
n th Coulomb branch, (|3.18|) , yields 

dE 40tt 2 m w 2 «-,« m w 2 ,„ 

—— « — = -75- « 1.616— ~- . (3.20) 

d£ v/2(r(l/4)) 4 i? 2 R 2 K J 

Hence, we see that the coulomb behavior of this branch persists all the way down to where 
it merges with the linear branch. 

We have thus seen that infinitesimal perturbations of the symmetry axis geodesic yield 
geodesies corresponding to a quark, anti-quark and n intermediate VF-particles experienc- 
ing a net attractive coulombic potential. We have also seen that there is a non-trivial 
threshold to go from n to n + 1 intermediate VF-particles. This means that we should 
view the symmetry axis geodesic as representing the process of adiabatically separating 
the quark and anti-quark. As the separation gets larger, the number of VF-particles that 
can be produced between the quarks increases. It is because of this increase in the W 
number that the energy does not undergo the Coulombic fall-off, but instead exhibits the 
confining behavior ( |3.11|) . 

A small perturbation away from the straight geodesic leads to a selection of a definite 
VF-number, and then further adiabatic separation of the quark and anti-quark preserves 
this number, but brings the state closer and closer to the BPS threshhold. It is easy 
to see that once a definite particle number has been selected, the confining behavior of 
E must at least turn over, and revert to a Coulombic decay. The reason is that if the 
energy increased without bound, then because the particle number is fixed, the energy will 
ultimately become far larger than the total mass of all the VF-particles between the quarks, 
and once again scale invariance will effectively be restored. We thus have the following 
picture of E(L) for all these states: the straight geodesic gives rise to and energy function 
E*{L) that is akin to the Cornell potential, starting as E*(L) ~ —1/L, and then crossing 
over to E*(L) ~ L. At finite, and ultimately regular intervals, the graph of E(L) versus 
L bifurcates: the confining potential, E*(L) continues, but a new branch, E n (L), drops 
away and corresponds to a quark, anti-quark and n VF-particles. This branch then follows 
the Coulombic behavior ( |3.18| ), shifted by the mass, ti^/tt, of n VF-particles. 
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3.3. QCD Configurations and Stability 



As we discussed in section 3.2, the physical picture for the (n + l) th coulomb branch 
is that of n VF-particles equally spaced between the quarks. This picture is confirmed by 
considering the energy density, p, along the string, which is given by ( |2.8| ). For a particle 
oscillating between the peaks, the QCD energy density gets very small as the geodesies 
approach the peaks, and becomes finite, or large as the geodesic passes across the symmetry 
axis. Indeed, as the particle moves higher and higher up the peaks it spends most of its 
time on the peaks, and very little time in between. Thus the QCD energy density becomes 
very tightly peaked, and if the particle executes n half-periods then there are n + 1 such 
peaks, the first and last being the quark and anti-quark, and the remaining n — 1 being 
VF-particles. For large n, the particle motion approaches the saddle and the oscillations 
become regularly spaced and so the VF-particles become very regularly strung out along the 
QCD string between the quark and anti-quark. Figure 4 shows plots of the string energy 
density for various values of L on the coulomb branch with six intermediate VF-particles. 
The plots show that the VF-particles quickly localize as the quarks are pulled apart. 
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Fig. 4: Plots of p vs. £ for the coulomb branch with six intermediate W's. p is 
in units of m\v 2 /R 2 and £ is in units of R 2 /mw- The quark separations are (a) 
L = 2.19R 2 /m w (b) L = 2.72R 2 /m w and (c) L = 3.37 R 2 /m w . The critical value 
for the appearance of this branch is L Pt 2.07 R 2 /mw ■ 



Thus these oscillating geodesies correspond to W-particles strung out like beads. For 
the geodesic that runs up to the saddle without oscillating, the VF-particles are smeared 
out along the length of the string, and as the oscillations increase the VK-particles localize 
to a greater and greater degree. Because of this attractive Coulomb potential between the 
beads on the string, we see that this configuration is very unstable, and so we may infer 
that the linear flux tube is also unstable. 
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We can also see this instability of the flux tube using perturbative analysis: Consider 
geodesies that run straight up the saddle and look for small perturbations 8x in the x 



direction. Using equations (|2.3| ) and (|2.6| ) we see that perturbation of the energy is 

2 dx< 



SE = - I t l( dSx ] 2 




(3.21) 



The variation of the length is 

tt = 2 / dA (^ fa+ ^) (3 ' 22) 

and so a quadratic variation in x can be compensated for by a linear variation in y in order 
that the length remain fixed. Thus, we see that the string configuration is unstable if the 
operator 

8 2 d 2 V 

C = ~ ~dx^ (3 - 23) 



has negative eigenvalue solutions. But we know from ( |3.12|) that a coulomb branch occurs 



precisely when this operator has a zero eigenvalue. Therefore, every time the linear branch 
passes by a coulomb branch, an eigenvalue of C changes sign. Hence if L n < L < L n+ i then 
the linear branch has n unstable modes. In terms of TU-particles, we expect an unstable 
mode for each particle and so the picture fits. 

To summarize, we see that the linear branch is unstable because the flux tube is 
energetic enough to pop out intermediate VF-particles between the quarks. The appearance 
of the VF-particles is marked by a localization of the energy density along the string between 
the quarks. The presence of these VF-particles then destabilizes the configuration because 
of a coulomb instability. We have seen in ( |3.19| ) and ( |3.2U| ) that the coulomb behavior of 
the n branch persists all the way down to L = L n , so we expect very little softening of 
this instability on the linear branch. 

When L < L±, then the branch of straight trajectories is stable and is in fact the 
only solution that extremizes the energy. In this case the Z2 symmetry is unbroken and 
the quark and anti-quark are in the fundamental of SU(N). When L > L\ then the 
stable branch is the first coulomb branch which is doubly degenerate. The degeneracy of 
this branch reflects the fact that the quark and anti-quark are fundamentals of one of the 
SU(N/2) groups. Thus, if we vary L between values that are less than L\ and greater 
than L\ such that we always remain on the stable branch, we see that there is a third order 
transition in that d 2 E/dL 2 is not continuous at L = L\. At this point the ^2 symmetry 
is broken. 
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4. Asymmetric and more general situations 

4-1. Trajectories off the Symmetry Axis 

Up to this point we have been assuming that the L>3-brane at infinity is on the 
symmetry axis defined by the two sets of D3-branes. For this configuration the heavy 
quarks are in the fundamental representation for one of the two SU (N/2) groups and both 
types of quarks have the same BPS mass. If the D3-brane at infinity is not on this axis, 
then the quarks will have unequal masses. If we denote the lighter quark by q and the 
heavier quark by q' ', then their mass difference is Am = <j> cos(0) / n , where 9 is the angle 
away from the D3 brane axis. So if 9 = degrees then the difference in quark mass is the 
mass of the W. 

We now consider trajectories that are at an angle 9 at asymptotic infinity. Quali- 
tatively, we expect to find behavior similar to the symmetric case for large L. That is 
we expect to find trajectories corresponding to the coulomb potentials for qq with some 
number n of intermediate TU-particles. As L is decreased, the energy of this configuration 
will match the potential energy for q'q with n — 1 VF-particles. Then as L is decreased 
even further this will match the potential energy for q'q' with n — 2 VF-particles. However, 
to make the next step to n — 3 VF-particles, we have to first jump down to qq with n — 2 
VF-particles, which releases a finite jump in energy because of the mass difference between 
the q and q' . Hence these last two geodesies will not be continuously connected. 

Hence, we expect to find a lone geodesic that corresponds to qq and which exists for 
any value of L. All other geodesies will come in bifurcating triple families. That is, the 
n th "triple" will correspond to qq with n + 2 VF-particles (n is even), q'q with n + 1 W- 
particles and q'q' and n W-particles. These three families of geodesies will be distinct for 
L greater than some critical amount. At very large L, the strength of the Coulombic force 
will be different for the different branches, and the energy curves will have an asymptotic 
separation of the VF-boson mass. As L decreases the energy curves, E(L), must merge, as 
outlined above, and E(L) must be singular (have a cusp) at the last merger. 

One can easily see how this picture of triples emerges from the symmetric picture 
considered in section 3. Consider the functions, i?*(L), and E n (L) described at the end of 
the last section. For n even the relevant geodesies ultimately run up one peak or the other, 
come to a halt and reverse their course, whereas for n odd, the geodesies never come to 
a halt and are symmetric between the peaks as in Figures 1 and 2. The energy functions 
E n (L), for n even, are thus doubly degenerate, corresponding to qq and n Ws, or q'q' 
and n W : s. This degeneracy is lifted when we split the degeneracy of the quark masses, 
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and the bifurcating curves E*(L) and E n (L) split into triples along the E n (L) curves for 
n even. 

Returning to the asymmetric situation, we can locate the critical values for the length 
at which mergers occur by examining the general behavior of the geodesies as a critical 
point is approached. As we saw in the previous section, critical points may appear if a 
geodesic becomes unstable due to a small oscillation. A singularity of this type is shown 
in Figure 5. In this figure we will denote the quark associated with the left-hand peak by 
q and the other quark by q' . We are considering geodesies with an incoming angle of 120 
degrees, and there is one in particular that rolls up the potential, stops, and then rolls 
back down. This geodesic is the limit of geodesies that roll up one peak, move to the other 
peak and then roll back down. The latter class of geodesic may be thought of as the state 
q'qW. The limit where this loop closes is a singular point. For this critical value of L, the 
q'qW trajectory merges into the qqWW trajectory. This second family of geodesies may 
be seen in upper part of Figure 6: these geodesies roll up the left peak, roll over to the 
right, and then reverse their course. In terms of QCD, the larger coulomb potential for 
qqWW compensates for the larger BPS mass of its constituents at the critical point. For 
L larger than this critical value one can be on either the qqWW or the q'qW family, while 
for L less than the critical value these two families are merged and are indistinguishable. 



Fig. 5: Geodesies with an incoming angle of 120 degrees that loop and cross the 
symmetry axis twice. For large L, the trajectory corresponds to a q'q pair and an 
intermediate W. There is a minimum value of L where the loop contracts to a line 
that retraces itself. At this value, the trajectory becomes identical with the trajectory 
for qq pair and two W -particles. 
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There is another type of singularity which is not really a singularity of the geodesies, 
but results in a singularity for the large N Super Yang-Mills. Suppose we have a family 
of geodesies parameterized by some quantity a. Then E and L are both functions of a. 
However, E and L can be multivalued and if so there is some point cr where d a E = and 
thus, by arguments given in section 2, d a L = 0. If we plot £asa function of L for these 
families of geodesies, we would find a cusp where the a derivatives vanish. Hence this is a 
singularity and in fact corresponds to the minimum L and E for these families of geodesies. 
Figure 6 shows a family of geodesies that have an incoming angle of 120 degrees and roll 
up the second peak, stop and then roll back down. As we discussed this family contains a 
transition from qqWW to a merger with q'qW. Continuing beyond this we can smoothly 
deform these geodesies so that the "rest point" of the particle moves down the peak and 
then moves back up it, finally representing a q'q' state. Clearly this family of geodesies 
will have a minimum value of L since L diverges as the rest point of the particle moves up 
the peak. At the minimum of L, the q'q' trajectory merges with the q'qW trajectory. Note 
that this critical point does not necessarily occur when the particle energy is minimized. 




Fig. 6: Geodesies with an incoming angle of 120 degrees that cross the symmetry 
axis twice. The trajectories stop at some point and roll back in the opposite direction. 
The figure shows two branches for a given L. For large L, one branch corresponds to 
qqWW while the other branch corresponds to a q'q' pair. The fifth trajectory from 
the left is the limit of the loop trajectories in Figure 5. The third trajectory from the 
right corresponds to the trajectory with minimum length. 

We can also see in Figure 6 that there is a trajectory in this class of geodesies where 
the particle energy /i is minimized. Near this point, ^ = l/y/2/j is roughly constant, so 
the potential is linear in L. Hence there is confining like behavior for trajectories in the 
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neighborhood of this point. These trajectories are the merged q'qW trajectories. This 
appears to be a general feature. That is, when a critical value L n is reached such that the 
qqW 2n trajectories merge with the q'qW 2 ™ -1 trajectories, then for L < L n the potential 
is roughly linear in L. Figure 7 shows the lowest coulomb branch and the first two triples. 
In the lowest triple, qqWW merges with q'qW which then merges with q'q' . Comparing 
Figure 7 with Figure 3, one sees that the effect of the unequal mass is to pull the degenerate 
branches in the symmetric case apart, leaving isolated sets of branches. 

2.5 

2 

1.5 

1 

0.5 



-0.5 

Fig. 7: Plot of E vs. L for the qq branch, and the two lowest sets of triples. As 
L is decreased the upper two branches in a triple merge merge. At even smaller L 
this branch merges with the bottom branch in the triple. To decrease L even further 
requires a finite jump in energy to a lower triple. 

Thus once one moves off the symmetry axis, the interesting features of the more 
symmetric situation survive, but in a more limited form. For example, there still appears 
to be a nearly linear potential for a range of separations, and in this range the number 
of Ws is changing, but the asymmetry limits the amount of "continuous" change in the 
particle number, and hence the behavior ultimately bifurcates into various Coulombic 
limits. 

The lone qq branch is obviously stable for any value of L since it is the lowest energy 
state. At L less than some critical value, L c , this is the only available branch. For L > L c 
a bifurcated triple branch appears. The lowest leg of the triple is the q'q' branch. While 
this branch has more energy than the qq branch because of the higher quark mass, this 
configuration is still a local minimum. To see this, note that in order for a q' to switch to 
a q, we need to pop a W near the quark, but the W will feel a net attractive force back 
to the quark. If we could pull the W far enough away from the quark, it would eventually 
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feel an attractive force toward the q' antiquark, and then these two would combine to form 
a q antiquark. Hence, we expect the q'q' branch to be stable, but the branch destabilizes 
when it merges with the q'qW branch, which occurs at the cusp at L = L c . There is then 
a finite change in energy as the configuration jumps down to the low branch. All other 
branches are unstable. 

It is now relatively easy to see the general structure of what will happen in breaking 
SU(N) to SU(Ni) x SU(N 2 ). First, there is still a saddle point, and thus there will be 
a "confining" geodesic running up to it. The confining geodesic will approach the two 
peaks of the potential V at some angle other than 90 degrees. One can view this angle 
as the asymmetry in the potential being compensated for by an asymmetry in the quark 
masses. There will be small oscillations about the confining geodesic, and much of our 
earlier discussion will go through. Indeed, it is easy to convince onself that the asymmetric 
potential, V, will result in families of geodesies that are deformations of the ones considered 
earlier. The corresponding conclusions about the quark-antiquark potential will thus be 
completely analogous. 

If one breaks SU(N) to p simple factors one can have a much richer saddle-point 
structure in the potential V. This will lead to different "confining geodesies," and to far 
more complex structure. In the 7L V symmetric breaking to (SU(N/p)) p , where the D3- 
branes sit on the vertices of a regular p-gon, one might expect that the Z p symmetric 
"confining geodesic" running along the symmetry axis of the p-gon would become more 
stable if the curvature in the saddle is softened. 

5. Final Comments 

N = 4 supersymmetric gauge theories are conformal and cannot confine, and yet we 
see geodesies that have some confining behavior. There is no real conflict here since the 
"confining geodesies" are classically and quantum mechanically very unstable, and the 
corresponding QCD state will destabilize by popping out VF-particles between the quarks, 
which are in turn swallowed by the quarks because of the coulomb attraction. In the end 
we are left with two quarks feeling a coulombic force. Of course linear flux tubes are 
not stable for real QCD either. If two quarks are pulled far enough apart, the flux tube 
destabilizes by popping a shower of hadrons out of the vacuum. 

On the more technical level, we have also found a simple relation for the interquark 
force in terms of a classical particle's energy. Through this we have shown that the geodesic 
prescription of |2||| always leads to an attractive force. It also gives a method for finding 
transitions in particle number, and for finding confining families of geodesies. 
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